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1. Introduction 

Circumstellar disks usually contain only a few percent of the total material 
going into a young star after the main collapse has stopped and the sur- 
rounding molecular cloud is cleared away. Yet the disks are of great interest 
to the study of star formation, perhaps as great as the stars themselves, 
because the disks may build planetary systems. The planet Earth contains 
less than one millionth of the mass of the Sun, but it is probably the most 
interesting body in the Solar System, certainly to us. Beckwith & Sargent 
(1996) argue that the currently known properties of disks are evidence that 
other planetary systems are common in the Galaxy and discuss the reasons 
for the interest in disk properties; that article provides a broad introduction 
to the subject. 

The purpose of this chapter is to provide a tutorial on how observations 
of the radiation from disks may be used to elicit their physical characteris- 
tics. In keeping with the spirit of the Crete meeting, the treatment is not a 
comprehensive review nor will it give a complete analysis of each method 
used to tease the disk properties from faint light observed with telescopes. 
Rather, the idea is to show that basic intuition about disk physics is eas- 
ily related to what is observed and provide a general foundation for the 
understanding of more elaborate theoretical calculations. 

Early disk models assumed that matter is confined to a very thin plane 
extending from the stellar surface to a sharp outer edge more than 100 AU 
from the star. The disk energy balance was attributed to accretion of matter 
through the disk. This oversimplified picture has been modified by a careful 
treatment of the underlying physics, and the more modern view is that the 
disk flares gently, often with an inner edge at some distance from the star, 
and is heated mainly by radiation as opposed to accretion. Some disks are 
surrounded by spheroidal "halos" that trap radiation and contain strong 
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outflows. Most of the young disks are accompanied by mass loss in columns 
along the polar axes that contribute to the total energy budget. Although 
it is not always possible to derive disk characteristics unambiguously from 
observations, most of the intuitive interpretations have been supported by 
increasingly better data and improved angular resolution images allowing 
us to separate the different components of a star/disk system directly. 
The article is organized along the following questions: 

1. What are the expected disk properties based on the theory of Solar 
System formation? 

2. How do we identify disks? 

3. How do we determine physical properties of disks from radiation? 

4. Do the observed properties show that disks are interesting? 

Because this article is a tutorial, some of the material is adopted from 
articles that I co-authored for Nature (Beckwith & Sargent 1996) and Pro- 
tostars and Planets IV (Beckwith, Henning, and Nakagawa 1999). 

2. The early Solar System 

It is generally agreed that a flat layer of gas and dust - a disk - orbited 
the early Sun and provided the material which later made up the Earth, 
Mars, Jupiter, and the other planets (Safronov 1969; Wood & Morfill 1988; 
Cameron 1988). The young Sun and the circumsolar disk were born from 
an extended cloud of gas and dust that was assembled from the detritus of 
dying stars and remnants of the early universe that collapsed under its own 
gravity. The material accumulated quickly onto the central proto-Sun but 
with enough residual angular momentum to prevent some from spiraling 
inwards - the exact proportion remaining in orbit is not known but should 
have been a considerable fraction of the total mass (Shu, Adams, & Lizano 
1987; Bodenheimer 1995). The average angular momentum of the collapsing 
region defined a rotation axis around which the orbits quickly stabilized, 
creating a disk with a thickness much smaller than its radius, at least within 
the regions now containing the giant planets. The formation of the stable 
disk probably occurred over about 10 5 years after the onset of free fall 
collapse (Shu et al. 1993), almost instantaneously in cosmic time. 

As the central proto-Sun evolved, the solid particles in the orbits settled 
to a dense layer in the mid-plane of the disk and began to stick together 
as they collided (Safronov 1969; Weidenschilling 1987; Mizuno et al. 1988). 
During the next 10 4 to 10 5 years, large rocks and small asteroids grew 
gradually from the small dust particles (Weidenschilling & Cuzzi 1993). 
When the gravitational pull of the largest asteroids was sufficient to attract 
neighboring pebbles and rocks, they grew even more rapidly to the size of 
small planets (Wetherill & Stuart 1993); gravity was important for bodies 
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more than 10 km across. The terrestrial planets are large accumulations of 
solid particles that grew from the collisions between these smaller bodies. In 
the outer parts of the disk, a few such solid cores became large enough (10 
Earth masses) to accrete gas (Mizuno 1980; Stevenson 1982), the dominant 
reservoir of mass, and gave rise to the giant gas planets (Wetherill 1990). 
Temperatures close to the proto-Sun were presumably too high to allow 
gas accretion. The planet building phase is thought to have taken between 
~ 10 7 and a few times 10 8 years, although the cores probably developed 
quickly, within the first 10 6 years or so. These timescales are not very 
well constrained by data, and it may well be observations of developing 
planetary systems around other stars that tell us how planets are really 
built. 

By analogy, we expect circumstellar disks to contain a few percent or 
more of the stellar mass, to extend at least 50 AU from the central star, to 
be relatively flat, and to be free from disruption for at least a few million 
years if they are to create the rocky cores needed to build the large planets. 
These characteristics certainly represent only a subset of the disks that ac- 
company star formation. Conceivably, a disk with substantially lower mass 
(10 -6 M Q ), size (a few AU), and lifetime (~ 10 6 yr) could create terrestrial- 
like planets suitable for life without the presence of gas giants. In principle, 
even larger, more massive disks could accompany the birth of stars much 
more massive than the Sun. Although we must keep an open mind about 
the characteristics the constitute a disk, the early Solar disk provides a 
framework to identify those disks that may become interesting for planet 
formation. 

3. How do we know that disks exist? 

Soon after the discovery that T Tauri stars - very young stars of approx- 
imately solar mass - had more radiation at infrared wavelengths than the 
photospheres should emit, Lynden-Bell and Pringle (1974) suggested that 
most of their peculiar characteristics might be explained by circumstellar 
disks. Their suggestion was based on the unusual spectral energy distribu- 
tions: the stars radiated too much ultraviolet light and too much infrared 
light at the same time. A disk could account for the ultraviolet light through 
emission from the boundary layer between the star and the inner edge of the 
disk, in which matter from the disk suddenly accreted onto the star, slow- 
ing down from Keplerian speeds to essentially zero speed so rapidly that 
the radiation temperatures are tens of thousands of Kelvin. The infrared 
light was radiation from the outer parts of the disk resulting from energy 
liberated as the matter slowly spiraled to smaller radii eventually to ac- 
crete through the boundary layer. One of their strong predictions was that 
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the long wavelength infrared radiation would follow a power law, F u oc V3, 
where F v is the flux density, and v is the frequency of the radiation. This 
result for the release of accretion energy through a disk is quite general. 

Even with no accretion, a disk will be heated by radiation from the star 
itself. The dust grains in the disk absorb stellar radiation and re-radiate 
in the infrared to maintain thermal balance. Remarkably, the spectral en- 
ergy distribution also follows a power law with the same exponent as for 
accretion over a broad range of wavelengths: F v oc between about 5 and 
100 fim depending on the luminosity of the star. At wavelengths short ward 
of about 3 /urn, the flux density stops increasing due to the inner radius of 
the disk, and at long wavelengths the power law becomes steeper due to 
the outer edge of the disks. The most complete treatment for a flat disk is 
given by Adams, Lada, and Shu (1988). 

When Lynden-Bell and Pringle made their suggestion, the long wave- 
length SEDs of disks could only be measured to about 10 fim. The SEDs 
could easily be explained over this limited spectral range by other distri- 
butions of dust near the stars. Although prescient, the predictions of the 
early disk theory went untested for nearly a decade. 

A few years after the SED calculations, Elsasser and Staude (1978) 
discovered that several young stars had rather high degrees of linear polar- 
ization in their optical light. They explained the polarization as scattered 
light from dust grains arranged symmetrically above and below a star that 
was obscured by a planar or toroidal distribution of dust oriented perpen- 
dicular to the line of sight and parallel to the direction of the polarization. 
These observations suggested that the dust distribution around the stars 
was axisymmetric and flattened relative to a spherical halo. 

The most striking demonstrations of axisymmetry in T Tauri stars are 
the well collimated jets seen in images of ionized lines. Mundt and Fried 
(1983) discovered the first of many young stellar jets in their image of 
HL Tau. The jets implied a strong axisymmetry near the stars, one com- 
ponent of which may be flattened disks perpendicular to the axis of the 
jets. In the same year, Cohen (1983) used the Kuiper Airborne observa- 
tory to measured the first far infrared SED of HL Tau, and he inferred the 
presence of a thin disk following the reasoning of Lynden-Bell and Pringle 
(1974). One year later, two groups observed HL Tau at angular resolutions 
several times better than previously possible and discovered elongated in- 
frared emission reminiscent of disks (Beckwith et al. 1984; Grasdalen et al. 
1984). Following the evidence accumulated for some time on this object, 
both groups suggested that the disk might be the progenitor of a plane- 
tary system similar to the early solar nebula. Subsequently, Beckwith et 
al. (1986) and Sargent and Beckwith (1987) made interferometric maps of 
the CO and 13 CO emission and showed that the gas, too, was elongated 
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as if in a disk; the direction of elongation was perpendicular to the jet and 
aligned with the major axis of the near infrared emission. The latter paper 
also observed the velocity field and interpreted it as gas orbiting the star. 
At this time, the evidence for a disk surrounding HL Tau appeared to be 
quite strong. 

These observations implied the distribution of material was disk-like but 
did not determine how much material was in a disk. The interferometric 
maps showed that HL Tau had continuum emission at 2.7 mm that appeared 
to be thermal emission from dust particles. With conservative assumptions 
about the radiative properties of the dust, the authors derived a mass for 
the HL Tau disk of ~ O.IMq with a considerable uncertainty. This mass 
is well above the minimum mass for the disk around the primitive solar 
nebula, 0.01 M , thereby strengthening the case for a protoplanetary disk 
surrounding HL Tau. 

At about the same time, astronomers analyzing the far infrared radia- 
tion from disks from the IRAS survey recognized that many T Tauri stars 
had strong far infrared excesses indicative of disks along the lines first noted 
by Cohen (1983). Rucinski (1995), Beall (1987) and Adams, Lada, and Shu 
(1987) recognized that the IRAS observations were most easily interpreted 
as disk emission, although there was some controversy about the applicabil- 
ity of Lynden-Bell and Pringle's flat, black disk calculations. Adams, Lada, 
and Shu (1988) showed that a class of these sources, termed "flat-spectrum 
sources" (see the next section) , required disks that were quite a bit warmer 
in their outer regions than predicted by the flat, black disk model. 

Two major advances followed in 1989 and 1990. Working with data from 
the IRAS sky survey, Cohen, Emerson, and Beichman (1989) and Strom 
et al. (1989) discovered that almost half the T Tauri stars in the Taurus- 
Auriga dark clouds had far infrared SEDs characteristic of heated dust, and 
the relatively low, visual extinction implied that the dust was in flattened 
distributions - most likely disks. This result meant that disks were common 
around young stars; HL Tau was not an anamoly. In the following year, 
Beckwith et al. (1990) found that a similar fraction of stars had emission 
at 1.3 mm, from which they placed an even stronger limit on the flattening, 
essentially showing that disks were the only viable explanations for the 
data. For the first time, they derived the distribution of masses for the 
disks around stars. Not only were half the T Tauri stars surrounded by 
disks, the great majority of the disks had sufficient mass to build planetary 
systems like our own. 

Most specialists in star formation had by then accepted disks as the best 
explanation for the accumulated data on T Tauri stars. But the evidence 
was largely indirect, and it was often possible to interpret the data sets 
with different theories for the distribution of dust near these stars. The 
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first pictures of disks taken with the Hubble Space Telescope demonstrated 
clearly that the dust distributions followed the theoretical pattern of a disk. 
The image of HH 30 by Burrows et al. (1996) shows a jet perpendicular to a 
dark disk seen edge-on with a gently flared shape, as expected for thermally 
supported disks (Kenyon & Hartmann 1987). Striking examples of disks are 
seen in silhouette against the Orion Nebula by O'Dell & Wen (1994) and 
McCaughrean & O'Dell (1996). Several more images of disks like HH 30 
have appeared recently (Stapelfeldt 1998; Padgett et al. 1999). Figure 1 
shows several images of disks taken with the Hubble Space Telescope. 
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Figure 1. These images of disks were observed with the Hubble Space Telescope. The 
upper three images are disks seen in silhouette against the Orion nebula (McCaughrean 
& O'Dell 1996). The bottom three images show disks close to edge on. The image of 
HH 30 (Burrows et al. 1996) shows flaring and illumination of cone-shaped cavities in 
the surrounding halo due to residual cloud material. The two images on the lower right 
are of main sequence disks: BD+31 643 (Kalas & Jewitt 1997) and (5 Pictoris (Kalas & 
Jewitt 1995, and references therein). 

The images leave no doubt about the distribution of matter: the parti- 
cles surrounding T Tauri stars are confined to flattened disks. Many other 
observations - very high resolution spectra (Hartmann et al. 1986), inter- 
ferometric maps of the gas (Dutrey, Guilloteau, & Simon 1994; Koerner & 
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Sargent 1995; Saito et al. 1995) and polarization maps (Bastien & Menard 
1990; Whitney & Hartmann 1992; Piirola, Scaltriti, & Coyne 1992) - are 
easily interpreted if disks are present and difficult to explain with other 
distributions of the matter. Disks exist. 

Compounded, the evidence for disks is overwhelming, but that does not 
mean that disks are present around every star or that individual stars may 
not have other distributions of matter. For those stars that do have disks, 
we turn our attention to their properties and how we observe them. 

4. Spectral Energy Distributions 

Disks contain a mixture of gas and dust. Early in a disk's evolution - the 
first few million years, say - there are enough small particles to absorb 
and reradiate light from the star. Dust warmed by the starlight radiates 
as a blackbody modified by the emissivity of the grains, which for small 
particles is approximately proportional to the frequency throughout the 
infrared region. The dust temperatures depend on the distance of a grain 
from the star. Because a continuous disk contains particles at a wide range 
of distances from close to the star's surface to several hundred AU or more, 
dust temperatures span several orders of magnitude. The result is a mixture 
of thermal emission across a broad spectrum from about 1 //m to more than 
100 fxm. 

The dust may be heated radiatively - by the central stars, for example 
- or by the gas via the loss of gravitational energy as it spirals through 
the disk to the star, called accretion heating. Accretion heating might be 
important in the earliest stages of disk formation and could dominate the 
total energy budget for a disk within the innermost 10 stellar radii or so, 
but it is probably unimportant for the infrared spectral energy distributions 
discussed in this section (Chiang & Goldreich 1996). For this reason, we 
will examine only radiative heating mechanisms for disks in this article. 

The study of accretion onto the star is a large topic in itself and well 
treated by other authors. The interested reader is referred to Bertout (1989). 
It has been suggested that disks may be heated by wave-driven accretion 
(Shu et al. 1990), where the waves are first moment modes (m = 1). This 
hypothesis has yet to find favor in the community, because it adds a com- 
plication to the physics for which there is no observational evidence, and 
radiative heating appears to be adequate to explain almost all observations. 
The subject is not closed, however, and it may be that dynamical heating 
is seen to be important when the angular resolution of telescopes becomes 
good enough to map the density distributions of disks. 

Young objects display a variety of different spectral energy distributions 
depending on their state of evolution. At the very earliest stages, the SEDs 
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are dominated not by the disks but by envelopes surrounding a central 
luminosity source most likely in free fall onto the star /disk system. After 
some time, these envelopes dissipate, either because all the matter is ac- 
creted onto the star or disk or because winds from the central objects blow 
away the remaining material. Lada (1987) invented a classification scheme 
to describe the sequence of SEDs that result from this process, a sequence 
generally used today to classify ages on the basis of SEDs. In the language 
of this scheme, we are concerned with Class II sources, sources in which 
the disk dominates the SEDs, and the envelope emission is negligible. The 
younger objects are treated by C. J. Lada in his chapter. However, Chiang 
& Goldreich (1998) show that some Class I sources can be explained as 
standard T Tauri star disks viewed at high inclination, and Padgett et al. 
(1999) recently published infrared images of Class I sources that show al- 
most all of them to be edge-on disks, so we must reserve judgment about 
the interpretation of individual objects based solely on their spectral class. 

4.1. FLAT, BLACK DISKS 

To understand the infrared spectral energy distribution of a disk, it is useful 
to calculate the simplest thin, black disk extending from the star 

to a distant boundary, specified by the outer radius, R m ax- Figure 2 is a 
schematic diagram of the star /disk system viewed exactly in the plane of 
the disk and showing the variables used for the following analysis. 




Figure 2. This sketch shows a star surrounded by a thin disk viewed in the plane of the 
disk. It extends from an inner radius, Rmin, to outer radius, Rmax (not shown), illumi- 
nated by a star of luminosity, L». The disk is axisymmetric in all physical parameters. 
A point at radius, r, in the disk is illuminated by a point on the star as shown with the 
variables labeled. The variable, x, is into the plane of the paper. 

The flux from the star that illuminates a point at radius, r, in the disk 
is given by integrating over the stellar surface as seen from the disk. It 
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is sufficient to consider only one surface (side) of the disk. If r is suffi- 
ciently large, so that the stellar diameter subtends a small angle, the flux 
is approximately: 



"ilium ~ -of 2sin0 / aT^dxdy, (1) 

r 2 Jo Jo 

2 A f R * fv\ rV R i-v 2 

K-Yrt, is) 



3 V r , 

where T* is the effective temperature of the star. The factor of sin 4> accounts 
for the illumination of the flat disk by a point on the star at distance, y, 
above the disk plane at distance, r, which is m (^) for large enough r. This 
flux is totally absorbed by a black disk. 

The disk radiates thermally as a blackbody, so: 

F rad = aT 4 (r), (4) 

from which we derive, 

1 3 

2\ 4 / r \ 4 



This result is nice, because it is easy to calculate the spectral energy dis- 
tribution with free parameters that are, in principle, observable. Assuming 
the star and disk are at distance, D, and oriented at angle, 9, to the line of 
sight (0 = corresponds to face-on viewing), the flux density is: 

1 f Rmax 

F„ = —p cos6B u [T(r)]2Trrdr (6) 

1 Anhu 3 ( kTo \ a 2 n f Xmax xdx 



D2 C 2 \ hv ) R ^ cos9 ; 3. . (7) 

u c \ in/ / Jx min exp{x<i) — 1 

where we define To = T(R m i n ), the other physical constants are labeled 
conventionally, and we substituted variables to create the definite integral. 
It is easy to evaluate (7) for any choice of inner and outer radii, but the 
beauty of this development is that for a wide range of wavelengths, the 
limits are essentially and oo, and the integral is a number. To see this 
result, note that the limits are: 

/ hv \i 



Xmm - KkT(R mir ^ (8) 
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For most wavelengths, the emission from the inner and outer parts of the 
disk is negligible; x m i n ps and x max « oo for all practical purposes. The 
definite integral is then: 



f 

Jo 



'o exp(x<i ) — 1 
and the spectral energy distribution becomes: 



« 2.576 (10) 
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Equation (5) shows how to determine T(R m i n ) for any value of R m in, if 
the stellar luminosity and source distance are known, and the disk orien- 
tation, 9, is observed with high enough angular resolution. The spectral 
energy distribution can be determined exactly. This predicted SED is the 
standard by which most SEDs have been measured to determine if a star 
is surrounded by a disk. 

The approximations made to arrive at (12) break down for high frequen- 
cies - short wavelengths - where the emission comes from near the inner 
disk boundary, R m in, arid at low frequencies where the emission comes from 
near the outer disk boundary, R m ax- There is a further modification that 
becomes important when the inner disk boundary is close to the stellar 
radius, R*. In this case, the star and disk mutually heat one another to 
increase the temperatures of both objects slightly - a kind of greenhouse 
effect - and the actual SED is more complicated: Adams, Lada, and Shu 
(1988) calculate this effect exactly for the interested reader. For our pur- 
poses, the approximations make very little error, since, as we will see, there 
are almost no examples of disks whose SEDs are fitted by this simple model. 

We derived (12) assuming a specific temperature distribution for the 
disk given by radiative heating. We can generalize the result for arbitrary 
power law temperature distributions for disks that have other ways of es- 
tablishing temperature distributions. Assume that T[r) = Tq{t j R m in)~ q , 
then (7) becomes: 



1 <iirhi/ 3 ( kT \? „ rxmax xc l x 

F v ~ 5— H— R minCOs9 



rxmax xc [x 

D 2 C 2 y hu J "rmn— exp{xq) _ V « ^ 

C 2 D- 2 coseT«R 2 mi y-l, (14) 
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where we used the approximation that x m i n ~ and x max m oo. The 
constant, C2 includes the definite integral evaluated numerically for values 
of q that allow convergence. 

Equation (14) is useful for SEDs that are power law in frequency. If the 
spectral index is defined as, F u oc u a , then (14) tells us that a = 3 — |. 
The spectral index, a is observed directly. Thus, it is always possible to 
derive the disk temperature distribution from the observed SED. This result 
generally applies to black disks, and they need not be very flat to produce 
this result. Even the absolute temperature at any radius may be derived if 
D and 9 are known. It is, therefore, possible to estimate disk temperatures 
with only modest uncertainties from the observed SED, assuming only the 
disk contributes (see Beckwith et al. 1990). 

The flat disk model is quite useful for developing intuition to make 
more complicated models easier to understand. At a wavelength, A, the 
thermal emission comes from a limited range of radii in the disk. The way 
to see this is to note that the integrand in (7) is maximum when x ~ 

0.5, meaning T(r) ps 1.7hc/k\, or r(A) « (0.6kT(R m i n )X/hc) 3 R m i n ~ 
4 

(A/A(i? m j n )) 3 Rmin- For a young star with = L Q , T* = 4000 K, and 
adopting R m in = R* ~ 2Rq, the wavelength reference is \{R m i n ) ~ 7 fim. 
At the outer edge of the disk, we can define a similar "reference" wavelength, 
that, for Rmax = 100 AU, say, is A(100 AU) = 7 mm. Emission from the disk 
dominates the SED throughout the far infrared, roughly 10 < A < 100 fim. 

It is also interesting to see what happens as the inner radius is increased. 
A disk with a hole in the center, such as seen around the main sequence 
stars Vega and (3 Pictoris, will not radiate at short wavelengths, because 
there is not enough material near the star. The effect is to lower the contri- 
bution from the disk in the near infrared without changing the far infrared 
radiation. SEDs characteristic of disks with inner holes are seen, and the 
changes brought about by a hole are pretty much the same for complicated 
disk geometries with flaring and radiative transfer as they are for the flat, 
black disk. Figure 3 shows a number of SEDs for different assumptions 
about the flat disks developed in this section. Keep in mind that some in- 
ner holes in the dust distribution will follow naturally from the melting of 
dust at ~ 1600 K near the stars. 

If a gap in the disk develops, it will also change the SED. Gaps are 
thought to be created when the first planets are assembled in disks. Un- 
fortunately, a gap will have only a modest effect on the SED, unless it is 
unusually large. The reason for this is that every wavelength gets contribu- 
tions from the disk over a fairly large range of radii, because of the width 
of the integrand in (7) and the slow variation of temperature as a function 
of radius. Therefore, gaps will be difficult to detect from the SEDs alone, 
unless they are very large. Figure 4 gives several examples. 
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Figure 3. These spectral energy distributions were calculated from the equations for a 
flat disk with different assumptions about the inner disk radii as denoted in the figure. 
The top half shows the separate contributions to the SEDs from the star and the disks 
with various inner radii. The bottom half shows the total SED consisting of the star 
alone, then the star plus the three disks with differing inner radii. 



The flat, black disk model makes a number of predictions that can be 
compared with observations: 

1. The infrared spectral energy distribution is a power law: F v oc us. 

2. The total disk luminosity, as measured by emission in excess of that 
from the stellar photosphere alone, is about -j that of the stellar lu- 
minosity. This fraction is exact for a disk extending from the stellar 
surface to infinity, and is in any case an excellent approximation for a 
disk extending a few AU from the star, since most of the luminosity 
comes from the inner radii. 



Circumstellar Disks 



13 



10 



vF v 



0.01 - 



o.i - 



1 - 




0.1 



1 



10 



100 



1000 



A,(um) 



Figure 4- These spectral energy distributions were calculated assuming gaps of different 
sizes are created in an otherwise continuous disk. The solid line is a continuous disk; the 
long dash is a disk with a gap between 0.75 and 1.25 AU; the short dash is a disk with 
a gap between 0.5 to 2.0 AU; the dotted line is a disk with a gap between 0.3 and 3 AU. 
Notice that the effect of a gap is slight until the gap cuts out almost a decade of radii 
from the disk. 

3. The spectrum of the disk should be smooth, if it is truly black. If there 
is some radiative transfer - i.e. if the disk has finite thickness - spectral 
features should appear in emission, because the disk atmosphere is 
heated from above. 

4. If there is additional luminosity very close to the star, as may occur if 
radiation is released by accretion of matter onto the stellar surface, the 
disk emission will also increase in response to this additional heating. 

It is easy to test these predictions. Figure 5 shows several observed 
spectral energy distributions of T Tauri stars with excess infrared emission 
compared with the predictions of the flat, black disk calculations. It is im- 
mediately obvious from the figure that the simple model does not fit the 
data. T Tauri SEDs generally have more emission at long wavelengths than 
the disk model permits. The infrared luminosities are, therefore, higher, 
and the slope of the long wavelength SED is flatter than the model. Mod- 
ifications of the model by increasing the inner radius, R m in, or decreasing 
the outer boundary, R m ax relative to the ideal model would exacerbate the 
differences with the data. The data show that if the T Tauri stars are sur- 
rounded by disks, the outer regions of the disks are hotter than they would 
be if heated by the central sources in a flat disk model. 

This problem was recognized as soon as the first SEDs were observed 
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Figure 5. Observed spectral energy distributions are plotted on calculations of SEDs 
using the flat, black disk model. In general, the actual SEDs have more excess infrared 
radiation than predicted by the flat disk model. 



(Rucinski 1985) and subsequently treated in an ad hoc manner by simply 
modifying the disk temperatures to empirically fit the SEDs (Adams, Lada, 
and Shu 1988). The ad hoc assumptions about disk temperature are un- 
satisfying in absence of a physical theory. The fiat, black disk model is not 
adequate to fit the great majority of data on circumstellar disks. 

The observations indicate that there is more far infrared radiation than 
expected from flat disks. Either the outer radii are hotter than would be 
the case for flat disks, or there are additional components - other reservoirs 
of dust - that can absorb and radiate along with the disk. To heat the outer 
radii, one must either increase the radiation field at large distances from the 
star or modify the disk in some way to absorb more energy from the star. 
There are models in the literature for all three approaches, additional dust 
components, increased radiation field, and increased absorption by the disk, 
and there is not yet a consensus on whether any or all of the possibilities 
are important in different objects. 
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4.2. MODIFYING THE DISK SHAPE: FLARING 

Disks should not be completely flat. Flaring will occur naturally in the 
outer parts of a disk due to the warming of the material. The disk should 
be in vertical hydrostatic equilibrium - i.e. normal to the surface - and ratio 
of thermal to gravitational energy increases with disk radius. In the outer 
parts, thermal energy fattens the disk where gravity is too weak to confine 
the material to a thin plane. For example, if the temperature as a function 
of radius, T(r) ~ r~ 3 ' 4 , as for a flat disk, and the vertical gravitational 
potential is dominated by the central star, E vert ~ — f ~ kT(r), then 
the scale height will vary with radius as, h sca i e ~ jj^j-r^/ 4 \ and the disk 
will flare - we ignore factors of order unity for this example. Figure 6 shows 
schematically what this means. 



Figure 6. The flaring of a disk occurs naturally for a disk in hydrostatic equilibrium. 
The disk mass is assumed to be negligible; gravity from the star acts to keep the material 
in a plane. The scale height of the disk increases with radius, because the thermal energy 
decreases more slowly than the vertical gravitational energy as radius increases. The 
vertical gravitational force, f vert , is shown as a component of the stellar gravitational 
force, f gravity The ray from the star shows the point at which short wavelength stellar 
radiation from the star is absorbed in the disk photosphere. The two other rays from this 
point show how the energy is reradiated into space and into the interior of the disk, thus 
heating the interior from the above. 

Flaring decreases the angle between the star and the normal to the sur- 
face, thus increasing the amount of light absorbed - the angle, </>,in (1) is 
increased. Kenyon and Hartmann (1987) recognized the importance of this 
effect as it applies to SEDs of disks. They noted that many SEDs can be 
fitted well by adopting some modest flaring of the disk surface and mod- 
eled the flaring necessary to produce good fits to extant data with black 
disks. Calvet et al. (1992) added the important effect of radiative trans- 
fer in a flared disk. They realized that the disk atmosphere would have a 
temperature inversion if heated by the star, and they explicitly calculated 
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models for flared disks with realistic atmospheres. Their numerical calcula- 
tions reproduced the main features of disk SEDs for a large number of disks 
and specifically predicted the presence of emission features, both molecular 
lines and resonant dust features, if the disks are irradiated from the central 
source. 

Chiang and Goldreich (1997) refined this idea by showing how the main 
physics, especially the flaring, could be understood analytically in a self- 
consistent calculation. They showed that realistic models of disk atmo- 
spheres do a good job of fitting data for most but not all disks. They also 
pointed out that because the dust particles that absorb the stellar radiation 
are in an optically thin part of the disk as far as re-emission is concerned 
and are small compared to the wavelengths of radiation, their equilibrium 
temperatures are warmer than they would be for purely black particles (cf 
Spitzer 1978). The upper disk atmosphere is warmer at larger distances 
from the star than predicted by the black disk calculations. Radiation at a 
specific wavelength must, therefore, have its origin over a larger area in the 
disk than expected from the black disk models. It will be easier to resolve 
these regions when very high angular resolution becomes available than 
predicted by some estimates (e.g. Bertout, Reipurth, h Malbet 1995). 

The results of these calculations produce two main changes to our think- 
ing about the physics of disks. First, the disk flares approximately as a 
power law with radius, where the height of the visible photosphere above 
the midplane, h p hot, as a function of distance from the star, r, is: h p hot ~ rJS 

4 

(Chiang & Goldreich 1997) or approximately ra. 

Second, in contrast to the black disk where radiation is effectively ab- 
sorbed and emitted from a black surface, light from the star at short wave- 
lengths - 0.5 to 2 fj,m, say - is absorbed in the visible photosphere at height, 
hphot, where it heats the dust. The dust then reradiates at longer wave- 
lengths, where it is optically thin above this layer and also downward into 
the disk until it is absorbed by dust in a lower layer that is opaque to the 
longer wavelength radiation. Because this longer wavelength radiation con- 
tains no more than half the energy originally absorbed in the photosphere, 
the inner layer is at a lower temperature than the photosphere. The result 
is a disk atmosphere in which the temperature increases vertically from 
the midplane outward, and most of the radiation seen by the observer is 
optically thin. Figure 6 is a shows the geometry of the radiative transfer 
and the different layers. 

The combination of optically thin and optically thick emission from the 
photosphere and interior, respectively, creates an SED that flattens out 
or rises toward long wavelengths. Figure 7 reproduces the calculation by 
Chiang and Goldreich (1997) showing how these models fit the data from 
the T Tauri star, GM Aur (see also Calvet et al. 1992). The flared disk model 
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can fit many of the disks in Figure 5 with modest adjustments to the free 
parameters. Since it is self-consistent and invokes no other components in 
the model except a star and a disk, it is a reasonable explanation for many 
observed SEDs. The superheated, optically thin layer obeys the scaling, 
F u ~ z/~ 5 / 7 (q = ||), and the interior, optically thick layer scales as, F v ~ 




Wavelength (urn) 



Figure 7. Figure 8 from Chiang and Goldreich (1997) showing how a flared disk with a 
photosphere reproduces one SED that differs substantially from a flat, black disk. They 
need a hole in the inner disk to account for the lack of disk emission short ward of about 
5 /J,m. 

The prediction of Calvet et al. (1992) and Chiang & Goldreich (1997) 
that irradiated disks produce emission features can be tested. A strong 
emission feature near 10 //m usually attributed to silicate resonance in dust 
was first observed for T Tauri star disks by Cohen and Whiteborn (1985) 
and more recently by Robberto et al. (1999). Figure 8 reproduces Figure 5 
from Robberto et al. (1999) showing spectra of 9 disks in the region of the 
silicate feature at 10 /iin, that can be compared with Figure 10 from Chiang 
and Goldreich (1997) showing the expected strength of emission features for 
a radiatively heated disk. The most beautiful example of emission features 
is given by Waelkens et al. (1998) and Malfait et al. (1998) from their ISO 
spectra of the Herbig Ae/Be star HD 100546 compared with the spectrum of 
comet Hale-Bopp. Figure 9 reproduces these spectra and demonstrates not 



18 STEVEN V. W. BECKWITH 




6 8 10 12 6 8 10 12 6 8 10 12 

A(/xm) A(yiim) X(/u,m) 



Figure 8. Observed spectra for the 9 stars in the sample from Robberto et al. (1999). 
The dotted line shows the continuum fitted to the data. The emission features are due 
to a resonance feature in silicate dust around 10 fim showing that these disks almost all 
have emission from optically thin dust, most likely in the superheated layers in the disk 
photospheres. 

only that emission features come from optically thin dust in the disk, but 
that these features indicate the dust composition is almost identical with 
that in the early solar system based on the comet spectrum. The results 
support the conclusion that disks are radiatively heated from above, as in 
the flared disk model. 

4.3. EXCEPTIONS: FLAT SPECTRUM SOURCES 

The disk SEDs that are most difficult to fit with the models described above 
are the flat spectrum sources (Adams, Lada, and Shu 1988). The princi- 
pal problem is that they emit more luminosity in the far infrared than at 
optical/near- infrared wavelengths, and so they defy explanations in which 
the disk temperatures are produced by radiation from the central star. Since 
a classical accretion disk produces an SED of identical shape to that of a 
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flat, black disk, T(r) ~ r _3 / 4 ) it does not help to increase the accretion 
rate. What is needed is to increase the temperatures in the outer disk, add 
another component for the far infrared radiation, or decrease the optical 
radiation in a very selective way; the optical/near- infrared colors are usu- 
ally inconsistent with pure extinction by normal interstellar dust. All ideas 
to understand the flat spectrum sources include additional components to 
the star and disk. 

A disk may be warmed by increasing the radiation flux on the outer 
parts if there is a way to scatter radiation onto the disk from the surround- 
ing region. Natta (1993) suggested that an optically thin spheroidal cloud, 
a "halo" , might scatter starlight in the outer regions of the disk to increase 
the energy absorbed and, therefore, the temperature. The spheroidal cloud 
may be a tenuous remnant of the molecular cloud from which the star and 
disk were created or it may be a transient part of the environment, such as 
residual infall or outflow in a wind. The halo acts like a greenhouse, trap- 
ping radiation at large distances from the star where the total scattering 
optical depth approaches unity. 

Calvet et al. (1994) expanded upon this theme to compute models of 
multi-component systems with a star, disk, outflowing wind, and a halo out 
of which the wind carves a cavity. These calculations can successfully repro- 
duce almost all observed SEDs in star /disk systems, but they include quite 
a few adjustable parameters, and they cannot be unique fits to the data. 
Calvet et al. (1994) argue that the flat spectrum sources are really closer to 
Class I sources in the current nomenclature than the Class II sources that 
concern us here, and so these models are perhaps less relevant to under- 
standing disk physics than the ones described in the previous subsection. 
Chiang & Goldreich (1998) fit Class I sources using only a disk viewed 
edge-on. They probably need an extra scattering component to account for 
the visible and near infrared light from some of these stars, but the number 
of assumption is, nevertheless, smaller. 

Calvet 's calculations demonstrate that realistic models of young stellar 
environments can fit the data in even the most truculent of SEDs. The data 
do not require new physics or surprising physical conditions near the stars. 
Figure 10 is an example of Calvet's model applied to the flat spectrum 
source, HL Tau. Her models require many assumptions about the distri- 
bution of matter near the stars, including infalling envelopes with cavities, 
winds, and non-spherical density distributions. In principle, most of these 
assumptions can be checked by observation, and it should be possible to 
determine how well the models correspond to reality. 

The introduction of additional parameters into these models does call 
into question our ability to derive physical parameters from the limited 
data sets provided by SEDs. Figure 5 shows that the number of indepen- 



20 STEVEN V. W. BECKWITH 



C/1995 01 (Hale Bopp) vs. HD 100546 




Wavelength(/im) 



Figure 9. This is Figure 5 from Malfait et al. (1998) showing the detection of emission 
features in the disk around the Herbig Ae/Be star, HD 100546, (solid line at top) as 
predicted if the disk is heated by radiation from the central source. It is compared with the 
spectrum of comet Hale-Bopp (dotted line underneath). Notice the close correspondence 
between emission features in the comet and in the disk spectrum, indicating that the 
particles in the disk are made from similar material as particles in the early solar system 
that made up the comet. 
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Figure 10. This is Figure 8 from Calvet et al. (1994) showing their fit to the flat spectrum 
source, HL Tau. Their model includes a disk, an infalling envelope with a cavity along 
the polar axes with an opening angle of lOdeg, and a scattering halo filling this cavity to 
produce the near infrared and optical light - shown as a separate line connecting open 
circles at the shortest wavelengths. This model successfully fits the SED at the expense 
of many adjustable parameters. 



dent points in a typical SED is limited. Often, there are at most five data 
points between 10 //m and 2 mm, the four IRAS measurements and a single 
millimeter-wave observation. The theoretical treatment of flat disks shows 
that each point derives its emission from a wide range of radii in the disk, 
and that rather large changes in the disk - large gaps, for example - impress 
modest changes on the detailed SED, details that may often be confused 
with noise when only a few points are available. The poor angular resolution 
of the IRAS observations means that far infrared data are susceptible to 
extended emission from the surrounding molecular cloud or a very tenuous 
but large halo surrounding the star. It is challenging, perhaps impossible 
to sort out the various contributions to the SED in this case. Therefore, 
the analysis of SEDs can give us only limited information about the disk 
upon which we can rely. The most important way to disentangle the vari- 
ous components will need very high angular resolution in the infrared (e.g. 
Padgett et al. 1999) and submillimeter. 



22 



STEVEN V. W. BECKWITH 



4.4. VERY LONG WAVELENGTH EMISSION: MEASURING DISK MASS 

At wavelengths longer than about 300 fim, these disks become optically- 
thin; you can see every particle in them. Therefore, the spectral energy 
distributions at these wavelengths are linearly proportional to the total 
amount of dust, and it is possible in principle to figure out the disk mass 
from the observations. The dust temperatures are almost always high enough 
to put this long wavelength emission on the Rayleigh-Jeans side of the 
Planck function, and the total emission becomes particularly simple to ex- 
press. 

The flux density, F u , from an optically thin disk at distance, D, is: 



Rmax 

B u [T{r)]T u (r)2irrdr, (15) 



where r u {r) is the optical depth. The optical depth can be written in terms 
of the surface density, X(r), and mass opacity coefficient, Kit as: T i /(/ 1 ) — 
K v T,(r). If the emission is in the Rayleigh-Jeans regime, B v m 2kTv 2 /c 2 . In 
this limit: 

2hl> 2 fRmax 

/ T(r)Z(r)2irrdr (16) 



■min 

2 



2k(T)u 2 „ r 

~ Kv c 2 D 2 M dust, (17) 

where (T), discussed explicitly in the next paragraph, is defined implicitly 
in (17). Regardless of the disk structure, the flux density is directly propor- 
tional to the particle mass opacity, k u , and the total particle mass, Md ust , 
through the integral. The absolute value of n v is difficult to determine accu- 
rately, although a good estimate may be made with reference to interstellar 
clouds. The frequency dependence of k u can be observed directly. 

If the temperature were constant, the integral is just the mass of the 
disk times the temperature: 

Z(r)27rrdr = M dust . (18) 



/ 

Jr.. 



The variable temperature has the effect of weighting the hotter (inner) parts 
of the disk more than the cooler (outer) parts. If the weighting function is 
not too steep, we can replace the integral in (17) with {T)Md ustl where 
(T) is some suitable average temperature that takes care of the weighting 
functions. We can calculate (T) explicitly for power law forms, T(r) = 
Toir/Rnin)-", and E(r) = T, Q {r / R min )-P-. 

I = / r S 2-Kvdr (19) 

J Rmin \RminJ 
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= 2vrr S i?Ln V :_ n _ n (20) 

z q p 

\ Rmin J 

« T M dust f§^)" 9 -^^, (22) 

V flmin / 2-q-p 

where we have assumed that 2 — g — p > for the last approximation and 

note that R m ax / Rmin ^> 1 > typi ca lly, Rmin ~ 

0.01 AU and R 

max ^ 

100 AU 

making this an excellent approximation. If 2 — q — p = or 2 — p = 0, one 
of the integrals is a logarithmic, and the approximation does not hold for 
2 — q — P < 0, but it is easy to work out the answers for those cases, too. 
For the case in (22): 

Thus, for q = 0.5 andp = 1, (T) « 2T{R max ), and F„ « K v ^^^ M dust . 

Using the analytical results from Chiang and Goldreich's analysis for the 
interior disk layer that produces almost all of the long wavelength radiation, 
we can work out (T) for a centrally irradiated, flared disk of the sort in their 
model. Noting that q = 6/14 and p = 3/2 in their calculations: 



( Rmax 
\ Rmin 



- 1 I 

(T) « T VJt "" n/ i -f (24) 

| Rmax J 2 ^ 14 

\ Rmin ) 

\tA^\\J^\ (25) 



I... „ //?„,„ V", /i? 



- ^T(i? max ) ^ In ^ (26) 
« 2.4T(i? maa; ). (27) 

For the numerical answer in (27), I assumed that (Rmax / Rmin) = 10 4 , but 
the very weak dependence of (T) on this ratio means that it is not a major 
uncertainty. This result is useful as a simple way to estimate disk mass 
from the temperature in the outer regions and the observed flux density. 
In most models including Chiang and Goldreich's, the temperature in the 
outer parts of the disk is around 20 K. This means that (T) rs 50 K for 
the purposes of estimating disk mass. In principle, it can be determined 
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directly from millimeter-wave maps of molecules in the outer disk. In any 
case, it is not very sensitive to the conditions in the disk for these models. 

The major uncertainty in mass estimates is in k u , not in the model pa- 
rameters. Theoretically, n v should be well determined for particles that are 
much smaller than the wavelength of observation, but it can vary substan- 
tially with changes in particle composition - silicates vs. iron, for example - 
and to some extent with particle size and shape. Attempts to determine k u 
from observations have given disparate results at millimeter wavelengths to 
date (Draine 1989; Beckwith & Sargent 1991), with variations of an order 
of magnitude published for ~ 1 mm wavelength. 

Theoretically, k u is expected to vary as a power of the frequency, u: 
K u = ^o( z/ / z/ o) /3 = fto(V'V)) - ^- A commonly adopted value at the time 
of this writing is 0.02(1.3 mm/A) cm 2 g _1 from Beckwith et al. (1990) and 
applies to total gas mass assuming Mg as /M(i ust = 100, but this value is by 
no means certain and could be wrong by a factor of five in either direction. 
Adopting this value, the disk mass can be written as a function of the 
continuum flux density at millimeter wavelengths as: 

F v ( D \ 2 f A \ 3 50 K 0.02 cm 2 g" 1 . . 

M disk = 0mM Q — — — —. 28 

Uy VIOOpc/ \1.3mmy (T) K Umin 



4.5. SUMMARY 

The spectral energy distributions of disks can reveal quite a bit of infor- 
mation about their physical properties. Some of this information is robust, 
and some is model dependent and must be regarded with skepticism. 

— Disks have broad thermal emission from the near infrared to the sub- 
millimeter, much broader than a Planck function (blackbody). The 
emission is often well fitted by a power law over an order of magnitude 
or more in wavelength. 

— The radial temperature dependence of the photospheric layer has an 
almost one-to-one relationship with the SED. A good estimate of the 
temperature as a function of radius, T(r), can be derived from the 
SED. 

— Typical disk temperatures are warmer than predicted by pure accretion 
heating or pure irradiation of a flat, black disk by a central source. 
Flaring of the disks is an economical means of explaining the warmer 
temperatures for many but not all disks. 

— Large gaps, including inner holes, will create "holes" in the SED, as 
well, although the effect of gaps is diminished because each wavelength 
derives contributions from a wide range of radii in the disk. 
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— Disks are expected to be optically thin (transparent) at millimeter 
wavelengths, and the total mass may be estimated from a measurement 
of the optical depth. 

— Disks with residual material in from their formation or from central 
outflows often display complex SEDs that require several components 
to reproduce. The SEDs alone are inadequate to demonstrate that the 
model distributions reflect the actual disk environment, but the model 
predictions can be checked from observations with high enough angular 
resolution. 

5. Properties of Disks 

It is not possible to image most disks in the manner of Figure 1. They are 
rarely situated in front of bright background sources, such as those in the 
top half of the figure, and they are normally not sufficiently edge-on to be 
seen in scattered light as they are in the bottom half of the figure. Detecting 
disks for the purposes of counting is normally done by looking for infrared 
radiation in excess of the stellar photosphere and assuming that any excess 
is due to disks. 

5.1. HOW MANY YOUNG STARS HAVE DISKS? 

As a general rule, it is easier to infer the presence of disks from radiation 
at wavelengths long ward of 10 fim, say, because that radiation is generated 
far from the star and requires more than trace amounts of dust to be 
observable. The first attempts to count disks in samples of young stars 
used far infrared and millimeter wave observations. Strom et al. (1989) and 
Cohen, Emerson, & Beichman (1989) found that nearly half of the young 
stars in their samples were detected in one or more of the IRAS bands - 
12, 25, 60, and 100 /im - and were likely to have disks. Beckwith et al. 
(1990) looked at similar samples at 1.3 mm detecting about 40% of them. 
These early studies asserted that disks were common among classical T 
Tauri stars. 

A variety of different surveys since then have generally supported the 
finding that disks are common, although the exact frequency depends on 
the sample and the selection criteria for young stars. It is much easier to 
detect light in the near infrared, at 2 /im, for example, where large arrays of 
sensitive detectors are readily available, but the amount of radiation from 
a disk is a much smaller fraction of the total - the star is bright, too - and 
very little dust is needed to produce an excess. Surveys of dark clouds at 
2 //m show disk fractions in anywhere from 10% to 90% of the sample. The 
high figure, 90%, is seen in the innermost regions of the Orion Trapezium 
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cluster (McCaughrean and Stauffer 1994), whereas the low figure includes 
large numbers of stars selected on the basis of x-ray emission (Walter et 
al. 1988). It is still controversial if x-ray emission is sufficient by itself to 
reveal very young stars, so the 10% figure should be regarded with caution. 
Recent surveys at 15 /mi with the ISO satellite give numbers in the 30 to 
50% range (Nordh et al. 1996; Nordh et al. 1998; Bontemps et al. 1998). It 
is reasonable to conclude that disks accompany young stars about half the 
time, and it may well be that all stars are born with disks. 

5.2. DISK LIFETIMES 

The lifetime of the dust is determined by counting disk fractions in young 
clusters with different ages. The first systematic attempts were made by 
Beckwith et al. (1990) at 1.3 mm and Skrutskie et al. (1990) at 10 /mi. The 
latter paper found that the sample fraction with disks decreased markedly 
between ages of a few million and ten million years. More recent studies 
using near-infrared excesses shows a steady decrease to almost zero at 10 
million years (E. Lada; this volume) indicating that the inner parts of the 
disks disappear over this time. A similar study at 25 and 60 /mi shows the 
outer parts of the disks disappear just as rapidly (Robberto et al. 1999; 
Meyer et al. 1999). 

Main sequence stars can also have tenuous disks, although the disk 
masses are many orders of magnitude less than their young counterparts 
(Backman &: Paresce 1993). These disks are presumably the remnants of 
those in the young stars, and they may be generated by the grinding of 
larger bodies - an asteroid belt, for example - orbiting the stars. Zuckerman 
and Becklin (1993) inferred a decrease of disk mass with time following 
approximately: M^isk ~ t~ 2 during the first 300 Myr. This result suggests 
a steady loss of small particles over almost the entire pre-main sequence 
lifetime of a star. 

Small particles should coagulate into larger ones on even faster timescales. 
The early solar nebula is thought to have created small planets within 
10 Myr. The observed disk lifetimes are, therefore, consistent with expec- 
tations for disks making planets, although there is no strong evidence that 
the particles are being accumulated into planet-sized bodies. 

5.3. DISK MASSES 

Most of the disk mass should be gaseous hydrogen and helium. The particles 
probably contain only 1% of the mass when a disk is first assembled from 
interstellar clouds. This ratio may change as the disk evolves. 

The particles are easy to observe, and the total particle mass can be 
estimated from the observed optical depth at millimeter wavelengths, as 
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we saw in a previous section. Millimeter wavelength observations of large 
samples show that typical disks are a few percent of a solar mass with a 
wide range. Figure 11 shows a histogram of disk masses for samples in the 
Taurus and Ophiuchus dark clouds (Beckwith et al. 1990; Andre et al. 1994; 
Osterloh and Beckwith 1995). These samples yield consistent results. New 
data (E. Lada 1999, this volume) indicate that the disks in the Orion star 
forming region may have lower masses on average, although the sensitivities 
of the observations are not yet good enough to make a robust conclusion. 
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Figure 11. The histograms show the distribution of disk masses among stars in the 
Taurus and Ophiuchus star forming regions determined by Beckwith et al. (1990)and 
Andre et al. (1994). 

Theories of the early solar system require between 0.01 and 0.1 M in 
the solar disk to create our planetary system. Disk masses around other 
stars fall in this range. It is clear that there is enough mass in a typical 
disk to create a planetary system like our own, but there is no proof that 
it will do so. 

If disks do make planets, there should be evidence of particle growth in 
the older disks. We know that particles are lost as the disks age. If the loss 
mechanism is accumulation into larger bodies, it seems likely that planetary 
systems are the final product of this evolution; if there is no particle growth, 
the disks may simply disperse into space or be accreted onto the central 
stars. 
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5.4. PARTICLE SIZES 

At submillimeter wavelengths and longer, k u is expected to vary as a power 
of the frequency or wavelength: n v oc oc \~@. For compact spherical 
particles smaller than the observing wavelength, (3 = 2 for metals and 
insulators under a wide range of conditions (Bohren and Huffman 1983; 
Emerson 1988). The power law exponent, (3, can be as small as 1 for cer- 
tain types of materials - amorphous carbonaceous material, for example - 
and even smaller over limited ranges, but it must go to 2 at long enough 
wavelengths to fulfill causality relations. It is expected that (3 ~ 2 for inter- 
stellar dust particles at long wavelengths. Careful attempts to measure both 
k u and (3 in interstellar clouds have, indeed, yielded [3^2 and values of 
Kv = (0.002- -0.004)(A/1. 3mm)" 2 cm 2 g" 1 (Hildebrand 1983; Draine and 
Lee 1984). If a mixture of particle types and sizes is present, the observed 
value is an average over the different constituents. 

Rocks, asteroids and planets are opaque to radiation at A ~ 1 mm, in 
which case (3 = 0. Pebbles with sizes of order 1 mm should have exponents 
with intermediate values: < (3 < 2. If dust grains grow large enough 
to put most of the mass in bodies larger than pebbles, they should have 
an observable signature at these wavelengths. Because the radiation does 
not penetrate far below the surface of a large body, meaning most of the 
material never interacts with it, the absolute value of the opacity per unit 
mass must decrease as the particles grow. 

Changes in the spectral energy distributions (SEDs) of disks relative 
to those typical of interstellar clouds are observed (Beckwith and Sargent 
1991; Mannings and Emerson 1994; Koerner et al. 1995). The flux densities, 
F u , of disks tend to fall more slowly at wavelengths longer than about 
400 ^m than those of interstellar clouds. The spectral index, a - where 
F v oc v a - is typically between 2 and 3, whereas a ~ 4 for the interstellar 
medium. Assuming optically thin emission in the Ray leigh- Jeans limit, the 
exponent, (3, is between and 1 for disks compared to the theoretical value 
of 2 for small particles in the long wavelength limit. The right hand part 
of Figure 12 shows the distribution of emissivity exponents, (3, determined 
from the observed spectral slope of disk emission near A ~ 1 mm. Many of 
the disks in this sample have (3 < 1. It is tempting to interpret this change 
as the result of particle growth. 

Figure 12 includes the distribution of /3's derived from model fits assum- 
ing the surface density is proportional to r~ 3 / 2 . Including the surface den- 
sity distribution tends to increase the derived value of the emissivity expo- 
nent relative to the observed spectral slope. For example at 1.3 mm, observa- 
tions of a disk with M D = 0.03 M , R D = 100 AU, and k u = 0.02 cm 2 g _1 , 
the average optical depth, f v = 0.24, and (3 « 1.36(q — 2). Any conclusions 
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Figure 12. The left hand part shows the spectral energy distribution of HL Tau from 
about 350 /im to 6 cm wavelength. The solid line fits the short wavelengths - thermal 
emission by dust - and the dashed line fits the long wavelengths - free-free emission 
from ionized gas. The long wavelength spectral index of the dust is 2.9, meaning (3 = 0.9 
with no corrections for optically thick emission (Wilner, Ho, & Rodriguez 1996). The 
right hand portion shows the distribution of /3's derived in two ways for a larger sample 
of disks. The "observed" /3's are derived directly from the spectral index as described 
in the text: (3 b s = a — 2. The "model" /3's are derived by assuming a surface density 
distribution as discussed in the text (Beckwith & Sargent 1991). 



about particle growth depend on knowledge of how much optically thick 
parts of the disks contribute to the spectral energy distribution at long 
wavelengths. It is probably for this reason that the samples of disks with 
measured a's at long wavelengths have not been enlarged much since the 
work of Beckwith & Sargent (1991) and Mannings k, Emerson (1994). 

To measure the optical depth, it is essential that the disk emission be 
resolved at wavelengths where it is relatively transparent so that the surface 
density distribution may be determined directly from the distribution of op- 
tical depth. Resolving the disks is time-consuming with the present suite of 
millimeter facilities. Several groups have now managed to resolve the emis- 
sion from one star, HL Tau, by carrying out novel experiments with existing 
millimeter- wave telescopes. Lay et al. (1994) and Lay, Carlstrom & Hills 
(1997) combined the CSO and JCMT into a single-baseline interferometer 
to observe at 0.65 and 0.87 mm and combined it with data from the Owens 
Valley millimeter interferometer at 1.3 mm. Mundy et al. (1996) extended 
the BIMA array to measure a size of the disk at 2.7 mm wavelength, and 
Wilner, Ho & Rodriguez (1996) used the VLA to get sub-arcsecond resolu- 
tion of the disk at 7 mm. Each group resolved the emission but with beam 
sizes of the same order as the disk size, and the sub-millimeter observa- 
tions had u — v plane coverage too poor to map the brightness distribution. 
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All groups measured the average brightness temperature and overall disk 
orientation and constrained subsequent model fitting much more than is 
possible from the SEDs alone. 

6. Are Disks Important? 

Disks are common around young stars, and they contain enough mass to 
make planetary systems like our own. At the time of this writing, sixteen 
planets orbiting stars other than the Sun have been inferred using radial 
velocity observations of the stars - although there is some controversy as 
to whether each of the objects is a planet according to usual definitions. 
Since we know some extrasolar planets exist and there are many disks, it is 
a small leap to suppose that disks make planets, as the theory asserts, and 
there are many planetary systems. This is an assumption, however, and 
proof is not likely to be forthcoming until we have actually observed some 
planets built up within the disks. 

If planetary systems are common and disks are the necessary precursors, 
disks are important for the understanding of planetary systems. They are 
likely to be our only means of filling in many missing links in our under- 
standing of early Solar System evolution. To understand our own origins 
and to assess the likelihood of other planetary systems, further study of 
disks around other stars presents one of the most fruitful avenues. 

On the other hand, the disks do not appear to be of great importance to 
the early evolution of the star after the main collapse has occurred. From 
the T Tauri stage until the main sequence, disks probably contain only a 
few percent of a stellar mass, and the amount of material accreted onto the 
star during this period is not likely to be sufficient to add significantly to a 
star's mass. It is possible that magnetic coupling between the disk and the 
star governs stellar rotation and determines the star's final rotation rate. 
There are some hints of this process but no proof. 

Many of us suspect that the very early stages of disk formation are 
crucial to the buildup of a star. During the initial collapse, there is almost 
certainly too much angular momentum in the molecular cloud to make a 
centrally condensed star unless a mechanism exists to shed angular mo- 
mentum during the collapse. A disk is a natural repository of such angular 
momentum. In this way, disks may be essential to a star's birth. But the 
disk angular momentum appears to be insufficient to take up that from the 
molecular cloud by the T Tauri phase (Beckwith et al. 1990), suggesting a 
large redistribution of angular momentum within the disk itself subsequent 
to the collapse phase. It is still difficult to study disks in great detail during 
and immediately after the collapse phase. There are not many candidates, 
yet - the entire phase is very short-lived - and the extinction toward the 
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very young disks is too high to allow even near infrared observations. Very 
high resolution millimeter wave images enabled by planned large interferom- 
eters (MMA/LSA) will likely change this situation and make observations 
of the physics of young disks routine. 
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